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We report that > 80% of the photons generated inside a photonic crystal slab resonator can be 
funneled within a small divergence angle of ±30° . The far-field radiation properties of a photonic 
crystal slab resonant mode are modified by tuning the cavity geometry and by placing a refiector 
below the cavity. The former method directly shapes the near-field distribution so as to achieve 
directional and linearly-polarized far-field patterns. The latter modification takes advantage of 
the interference effect between the original waves and the reflected waves to enhance the energy- 
directionality. We find that, regardless of the slab thickness, the optimum distance between the 
slab and the reflector closely equals one wavelength of the resonance under consideration. We have 
also discussed an efficient far-field simulation algorithm based on the finite-difference time-domain 
method and the near- to far-field transformation. 

PACS numbers: 42.70.Qs, 41.20.Jb, 42.60.Da 
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I. INTRODUCTION 



Multitudes of optical resonant modes can be created 
by introducing a structural defect into a perfect pho- 
tonic crystal (PhC) having photonic band gap (PBG). 
Many researchers have been interested in designing novel 
nanocavities based on the PhC in the hope of controUing 
the flow of light in unprecedented ways.|l| Especially two- 
dimensional (2-D) PhC slab structures have been widely 
studied by utilizing well-established standard fabrication 
processes. In the horizontal directions of such 2-D PhC 
slab structures, photons can be spatially localized by 2- 
D photonic bandgap effects. Additionally in the vertical 
direction, photons are confined rather efficiently through 
total internal reflection. [3, Q Air-hole triangular lattices 
are widely adapted as a basic platform of the 2-D PhC 
because of their large photonic bandgap. [2| 

After the first demonstration of the PhC nanocav- 
ity laser based on a single missing air-hole by Cal- 
tech group, Q various cavity structures have been 
demonstrated, showing rich characteristics such as 
nondegeneracy,|3 low-threshold, high Q/V value, |3,ll| 
and so on. Recently, 2-D PhC cavities have drawn 
much attention as a promising candidate for cavity 
quantum electrodynamics (CQED) experiments report- 
ing vacuum Rabi splitting 9] and high-efficiency single 
photon sources lOj. For CQED experiments, various PhC 
cavities with fine structural tunings have been proposed 
to achieve good coupling between the cavity field and 
an emitter placed inside the cavity. 11, 12] Park et al. 
demonstrated the first electrically-driven PhC laser based 
on the monopole mode that has an electric-field null at 
the cavity center where the current-flowing posttuj is 
placed. The 2-D PhC structure can also be used for in- 
plane integrated optical circuits where PhC cavities and 
PhC waveguides are coupled[ij to supply added func- 
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tionalities. By using a PhC-based channel add-drop fil- 
ter, one can drop (or add) photons horizontally (or verti- 
cally) through the cavity. [1^ [la, [l3 In this way, various 
functional PhC devices have been designed. 

In this article, we shall discuss useful design rules 
for the vertical out-coupling of PhC light emitters. To 
achieve efficient out-coupling into a single mode fiber 
using conventional optics, the far-field of the PhC res- 
onator should be both well-directed and linearly polar- 
ized. Gaussian-like far-field emission profile is generally 
preferable to have good modal overlap with the funda- 
mental mode of a fiber. Moreover, one needs to find 
ways to overcome the strong diffractive tendency of a 
wavelength-small cavity. We begin by looking for pos- 
sible remedies through the defect engineering approach 
(tuning the structure near the defect). Then, we will 
investigate effects of a bottom reflector on the far-field 
characteristics. The surface of high-index (~ 3.5) sub- 
strate is regarded as a bottom reflector at flrst. Note 
that such a reflector is practical since it can be natu- 
rally embedded in most of the wafer structures, [l^ [13 
One can also think of a highly-refiective Bragg mirror 
made of alternatin g la yers of high- and low- refractive in- 
dex dielectric slabs[23 as an alternative bottom refiector. 
We will investigate both the defect engineering approach 
and the bottom reflector effects to engineer the far-field 
radiation characteristics. A simple and an efficient far- 
field simulation method based on the 3-D finite-difference 
time-domain (FDTD) method and the fast Fourier trans- 
form (FFT)[2l|,|2i| wiU also be discussed. 

This article is organized as follows. In Sec. |nl vari- 
ous characteristics of 2-D PhC slab resonant modes are 
reviewed. In Sec. IIIII we will describe the far-field sim- 
ulation method and show simulation results. In Sec. II VI 
we will discuss how the defect engineering method can be 
applied to the PhC resonant mode to realize directional 
emission. Then, effects of the presence of a bottom re- 
flector will be discussed. Finally, we discuss how the two 
proposed schemes can be combined to give both direc- 
tional and linearly-polarized emission. 



II. RESONANT MODES IN A MODIFIED 
SINGLE DEFECT CAVITY 



In this section, we review some important characteris- 
tics of a few possible resonant modes in a modified single 
defect PhC cavity. 18] The 3-D FDTD method was em- 
ployed to analyze optical characteristics such as the mode 
profile, the Q factor, the resonant frequency, and so on. 
A typical structure of the modified single defect cavity 
and various resonant modes are shown in Fig. ^^a). A 
free-standing PhC slab structure is assumed as a plat- 
form for a triangular PhC lattice. The refractive in- 
dex of the slab is chosen to be 3.4, which corresponds 
to that of GaAs around 1 /im. For real applications, 
one can insert active semiconductor layers, such as the 
In As quantum dots, in the middle of the slab. The 
air-hole radius and the slab thickness are chosen to be 
0.35a and 0.50a, respectively, where a represents the lat- 
tice constant of the triangular PhC. The choice of the 
slab thickness is important to support only fundamental 
modes of the slab (TE-like modes) . 23] Without modifi- 
cation on the nearest air-holes around the defect, only 
the doubly-degenerate dipole mode can be excited. |j,|2J] 
However, by reducing the radius of the six nearest air- 
holes as indicated in Fig. ^a), other resonant modes 
within the PEG can be excited. This type of modifi- 
cation was first proposed by Park et al. so as to find a 
nondegenerate monopole mode used for a photo-pumped 
low-threshold PhC laser. (^ Especially, as indicated by 
Ryu et al., the hexapole mode can have a very high Q/V 
value due to the whispering gallery nature of the sym- 
metric field distribution. 8] For those modes having an 
electric-field intensity node at the cavity center, intrinsic 
properties of those resonant modes is not changed appre- 
ciably even after the introduction of a small post at the 
cavity center. This was a crucial design rule applied for 
the electrically-pumped PhC laser. [l^ One disadvantage 
of these resonant modes is that the vertical emission is 
inherently prohibited due to the destructive interference 
of the anti-symmetric field patterns. 

Resonant frequencies and Q factors of possible reso- 
nant modes in the modified single defect cavity are sum- 
marized in Table ^ Generally, the total radiated power 
(l/Qtot) can be decomposed into a vertical contribution 
(l/Qvert) and an in-plane contribution (l/Q/iorz)-|2J] 
One usually refers to 1/Qvert as the inherent optical loss 
of the resonant mode, because the in-plane loss {1/Qhorz) 
can be arbitrary reduced by increasing the number of 
PhC layer surrounding the cavity. In fact, the Q fac- 
tors listed in Tableware Qtot, which have been obtained 
by using a sufficiently large horizontal computational do- 
main (I6a X I6a). The optical loss of the resonant mode 
is closely related to the in-plane field distributions. For 
graphical illustration of the vertical loss, the momen- 
tum space intensity distribution'25| obtained by Fourier 
transforming the in-plane field {iFTiE^)]"^ + iFTiEy)]"^) 
is calculated. From the in-plane momentum conserva- 
tion rule, only plane- wave components inside a light-cone 
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FIG. 1: (a) Structure of a modified single-defect cavity. Here, 
the six nearest neighbor holes are pushed away from the cavity 
center after reducing their radii from r = 0.35a to rm = 
0.25a. (b) Various resonant modes in a modified single-defect 
cavity. Dipoles and quadrupoles are doubly degenerate while 
a monopole and a hexapole are nondegenerate. We plot Hz 
components of the each resonant mode, where the electric- 
fields are polarized in the x-y plane (TE-like modes). 



TABLE I: Normalized frequencies {un = a/A), Q factors, and 
theta-polarized percentages of the resonant modes shown in 
Fig-Hb). 



Pe/Ptotal{%) 



Dipole 

Quadrupole 

Monopole 

Hexapole 



0.2886 
0.3179 
0.3442 
0.3149 



14 900 
45 000 
11 000 
168 000 



46 
75.9 

1.2 
74.5 



[fc^ + ky — (oj/c)'^] can couple with propagation modes. 

As can be seen from the momentum space intensity 
distributions of the resonant modes in Fig. |21 the major 
fraction of the resonant photons live outside the light- 
cone, which implies an efficient index confinement mech- 
anism. The hexapole mode has negligible plane-wave 
components inside the light-cone, hence a high Q fac- 
tor of ~ 168 000. Among four possible resonant modes, 
the monopole mode and the dipole mode have relatively 
low Q factors {^ 10 000). Due to the presence of dc com- 
ponents (fca; — ky — 0), the dipole mode shows relatively 
large radiation losses. Although the controlled vertical 
emission is essential for good out-coupling, one should 
not give up the quality factor too much. In fact, several 
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FIG. 2: MomentuiTi space intensity distribution {\FT{Ex)\^ + 
|-FT(_Ej,)|^) of the resonant modes in a modified single-defect 
cavity. The dotted white circle in each graph represents a 
light-line defined by k'^ + ky = (uj/c)^. The major fraction 
of the plane-wave components is outside the light-line, which 
implies the efficient index confinement mechanism. 



groups have reported high Q/V PhC cavities with direc- 
tional emission. jTllllL Il4l| Later, we will investigate how 
one can perform additional fine tunings on the modified 
single defect cavity to obtain both directionality and high 
Q/V. 

Finally, we would like to mention how the structural 
symmetry of the PhC cavity can affect far-field char- 
acteristics. The modified single defect PhC cavity de- 
picted in Fig. n retains complete six- fold symmetry. The 
group theory tells us one can classify all the resonant 
modes into nondegenerate modes and doubly degener- 
ate modes depending on the rotational symmetry; the 
monopole mode and the hexapole mode are nondegener- 
ate and the dipole mode and the quadrupole mode are 
doubly degenerate. |2a l27l | Fig. El shows E^ and Ey field 
distributions of the resonant mode. First, let us look at 
the case of one of the degenerate dipole modes. In E^ 
field distribution, due to the odd reflection symmetries 
by both a; = and y = 0, the E^ vertical emission is 
prohibited. However, in Ey distribution, since all the 
reflection symmetries are even, there remains a dc com- 
ponent which contributes y-polarizcd vertical emission in 
the far-fleld. In the cases of the monopole mode and the 
quadrupole mode, at least one of the reflection symme- 
tries are odd, thus the vertical emission is always pro- 
hibited. In the case of the hexapole mode, E^ field dis- 
tribution shows even reflection symmetries by both axes 
(x = and y = 0). Considering its perfect six- fold sym- 



metry. Ex field components must be completely balanced 
to give a zero dc component. Thus, E^ field distribution 
reveals that such delicate balance can be easily broken 
by introducing a small structural perturbation, resulting 
in nontrivial vertical emission. We will discuss this type 
of defect engineering in Sec. llVI 



III. FAR-FIELD SIMULATION 
A. Far-field simulation method 

To obtain the far-field radiation pattern, we should 
deal with the radiation vectors that are sufficiently far 
(r ^ A) from a light emitter. [2 Sj Direct application of 
the FDTD method would be difficult because of the lim- 
ited computer memory size and time. Here, we shall 
explain how one can efficiently obtain the far-field radia- 
tion by combining the 3-D FDTD method and the near- 
to far-field transformation formulae. This method has 
been well known in the field of electronics since 1980s 
dealing with antenna radiation problems. Vukovic et al. 
applied the FFT-based far-field simulation algorithm to 
their Q factor optimization of the PhC cavity mode. |23| 

We explain the basic concept of the efficient far-field 
computation in Fig. 4. Here, we divide the whole cal- 
culation domain by a horizontal plane located just above 
the PC slab. The upper domain is bounded by an infi- 
nite hemisphere. All field components (E and H ) are 
assumed to fall off as 1/r, typical of radiation flelds. It 
is also assumed that in-plane field components (E^, Ey, 
Hx, and Hy) detected at the horizontal plane decay ex- 
ponentially as the horizontal distance p. According to 
the surface equivalence theorem, for the calculation of 
the fields inside an imaginary closed space fl, the equiv- 
alent electric (J^) and magnetic (M^) currents on the 
surrounding surface can substitute all the information on 
the fields out of the fl, where 3s and Mg are calculated 
by using the following formulae. 



Js = n X H = —xHy + ijHx 
Ms = — n X E = xEy — yEx 



(1) 
(2) 



Here, n is a unit normal vector on the surface. Neglect- 
ing the fields on the hemisphere, only the in-plane field 
data detected at the horizontal plane will be used for the 
far-field computation. Then, these equivalent currents 
are used to obtain the following retarded vector poten- 
tials. 



A = /io / — : dS 



F = er, 
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(3) 

(4) 



In the far-field regime, {r Ki R — pcos^), the above 
formulae can be simplified as 
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FIG. 3: Ex and Ey field distributions of the resonant modes in a modified single-defect cavity. For the doubly degenerate 
modes such as the dipole and the quadrupole, we plot one of the degenerate pair (the left ones among degenerate pairs shown 
in Fig. p. 
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FIG. 4: Geometry used for the far-field computation. The 
whole calculation domain is divided by the horizontal plane 
just above the photonic crystal slab. According to the surface 
equivalence theorem, the near-field components {E^, Ey, Hx, 
and Hy) at the horizontal plane give sufficient information on 
the far-field pattern in the upper domain (fl). 



^ — ikR p p — ikR 

^ — ikR. p p — ikR 

F = e„— — / M,e^^-''™^*d5 = £«— ^L (6) 

Using a relation fcpcos^ = kxX + kyy, the radiation 
vectors N and L are simply related by 2-D FTs. 



L= I M.se'^P""'''^ dS = FT{M.s 



(7) 
(8) 



X and y components of the radiation vectors can be 
represented as 
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(12) 



Then, far-field radiation patterns are obtained by cal- 
culating time-averaged Poynting energy per unit solid 
angle. |28i] 
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(13) 
where rj is the impedance of free space. 

One can easily implement this algorithm into the 
FDTD code. During the FDTD time stepping, the in- 
plane field data (near-field data) that contain complete 
information on the far-field are obtained. Once the 
FDTD computation is completed, the radiation vectors 
(N and L) are calculated by the FTs. Then, these ra- 
diation vectors are used to obtain far-field patterns ac- 
cording to Eq. p3l) . One can use the 2-D FFT algorithm 
for efficient FT computations [see Eqs. (l^- (|12|l ]. In our 
FDTD simulation, the size of the computational grid was 
a/20 (a is the lattice constant), and the size of the hori- 
zontal calculation domain were 16a x 16a. Thus, each in- 
plane field data will contain approximately 320 x 320 data 
points. In the FFT algorithm, N x N input data I{x, y) 
are transformed to the output matrix data 0{fx, fy) with 
the same N x N size, where / is fc/27r. 29] Assume that 
N unit cells in the I{x, y) data correspond to a spatial 
length of L wavelengths (A). The frequency step A/ 
is 1/(LA) = fo/L and the maximum frequency is thus 
Nfo/L. Thus, the value f — fo oi the light-cone corre- 
sponds to the Lth point in the frequency space. Note that 



only the wavevector components lying inside the light- 
cone whose radius is L cells can contribute to the outside 
radiation. To increase the resolution of the far-field sim- 
ulation, one should increase the horizontal size of the 
FDTD domain. Fortunately, there is an easy way to do 
so without performing an actual calculation in the larger 
FDTD domain. We can enlarge the input matrix data by 
filling null (0) values, for example, into 2048 x 2048 size 
matrix. Such arbitrary expansion is legitimate since most 
of the near-field energy of a PhC cavity mode is strongly 
confined around the central defect region. When the nor- 
malized frequency of a PhC cavity mode is a/A ~ 0.35 
(typical of PhC defect modes), L = 2048 cells = 102.4a 
~ 36(A). Thus, the resultant far- field data resolution is 
A/ - hIL - /o/36. 

Finally, we would like to mention one more before pre- 
senting calculated far-field patterns. All E and H data 
described in the above formulae are phasor quantities. 
Thus, when we detect the in-plane field data, we have to 
extract both amplitude and phase information, because, 
in general, phases will change differently at each position. 
The situation containing the bottom refiector is the case, 
where the reflected waves can no longer be described as 
a simple standing wave of E(r, i) — |E(r)|e"^*. Assuming 
a single mode uj, all field vectors in the simulation should 
have the following form. 



E(r, t) ~ E(r) exp(ici;t) — |E(r)| exp(i</7(r)) exp(ici;t) 

(14) 

An efficient method to obtain the exact phasor quanti- 
ties E(r) is to incorporate recursive discrete FTs (DFTs) 
concurrently with the FDTD time stepping. 
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FIG. 5: Calculated far-field patterns of the resonant modes 
depicted in Fig. 3. All the far-field data {x, y) are represented 
by using a simple mapping defined by a:: = 6' cos 4> and y — 
Ssin</) (The radius of the plot corresponds to 6*.). Along with 
total intensity (|Ep) patterns, we also display theta-polarized 
intensity (|_Bep) patterns and phi-polarized intensity (I^E^p) 
Ex\i,j ■ cos{uj ■ nAt) (15) patterns. Intensities of \Ee\'^ and \E^\'^ are normalized by the 
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sin(w • nAt) (16) 



maximum value of both data. 



Here, to obtain real (Exicai) and imaginary (Exi,„ag) com- 
ponents of Ex field, simple recursive summations are per- 
formed on each grid point {i,j) at the FDTD time n. 
Note that the above formulae correspond to the follow- 
ing Fourier transformation. 



Ex{r,iu)^ / Ex{v,t)e'^'dt 



(17) 



In general, this summation should continue until the 
field values 'ring down' to zero. However, in the case 
of high Q factor resonant mode, by carefully exciting 
a dipole source (source position and frequency), small 
FDTD time steps (usually 20 000 ~ 30 000 FDTD time 
steps, or equivalently ^ 200 optical periods) are sufficient 
to obtain nearly exact phasor quantities. In the next 
section, we will show calculated far-field patterns. 



B. Far-field simulation results 

To calculate the far-field pattern, we have used the 3-D 
FDTD method and the near- to far-field transformation 
algorithm which has been outlined in the previous subsec- 
tion. The horizontal calculation domain of 16a x 16a and 
the spatial grid resolution of a/20 were used. Calculated 
far-field patterns are shown in Fig. [S] All the far-field 
data {x, y) shown in this article are represented by using 
a simple mapping defined hy x = 9 cos cf) and y = sin (f>. 
Depending on the rotational symmetry that the each res- 
onant mode has, the far-field pattern also shows several 
regularly spaced lobes, which in turn explain the labeling 
of these modes; Number of lobes in the far-field pattern 
are four (six) in the case of the quadrupole mode (the 
hexapole mode) . Ignoring weak intensity modulation laid 
upon dominating circular emission pattern, the monopole 
mode shows two concentric ring-shape patterns, which 



show indeed a monopole-like oscillation. In fact, emission 
patterns of the hexapole mode and the monopole mode 
does not bear perfect six- fold symmetry, e.g. two hori- 
zontal lobes look bigger than the others in the case of the 
hexapole mode. These discrepancies are due to rectan- 
gular grids of the FDTD computation cell which are not 
well matched with the six-fold symmetric case. (Later, 
we will explain more in detail.) Along with total intensity 
(|Ep) patterns, we also display theta-polarized intensity 
(|£'ep) patterns and phi-polarized intensity QE^l"^) pat- 
terns. Note that such polarization-resolved far-field data 
can be easily obtained in the real measurement setup 
(a solid angle scanner) by placing a polarizer in front 
of a photo-detector which scans the whole hemispherical 
surface. |3(TJ 

First, let us look at the results of the quadrupole mode. 
Comparing with the total intensity power, it can be de- 
duced that most of the emitted power is theta-polarized. 
Also in the case of the hexapole mode, theta-polarized 
emission seems to dominate phi-polarized one. Table |2 
summarize percentages of the theta-polarized power and 
the phi-polarized power. As expected, the quadrupole 
mode and the hexapole mode are almost theta-polarized 
(> 70%). However, the monopole mode shows opposite 
result; phi-polarized emission dominates by ~ 99%. As 
indicated by Povinelli et at, one may create an artifi- 
cial magnetic emitter by using the monopole-like mode 
which shows the same emission characteristics reminis- 
cent of those of the magnetic multipole.|3H Consider- 
ing that electric (magnetic) dipole radiation is theta- 
polarized (phi-polarized), [28J such polarization charac- 
teristics seem to be closely related with composition 
of the electric multipoles and the magnetic multipoles. 
Thus, the so called multipole decomposition could pro- 
vide other perspectives for various radiation properties 
of PhC nanocavities along with the plane-wave decom- 
position (FT based analyses). 

From polarization-resolved far-field patterns, it looks 
as if the dipole mode has a certain singularity near the 
north pole {9 = 0). This feature comes from the strongly 
2/-polarized dipole mode. Upon inspection of the region 
around the vertical, one can observe that directions of the 
electric field vectors indeed predominantly lie in the y di- 
rection. Fig.l^a) shows the a;-polarized emission pattern 
( I £a; P ) and the y-polarized emission pattern ( | i<Jy p ) . The 
result shows the domination of the y-polarized emission 
over the x-polarized one. _32j Thus, except the inherent 
degeneracy, the dipole mode itself can be a good candi- 
date for the vertical emitter. 

Finally, we would like to mention the origin of the 
unevenness for the six lobes in the far-field patterns 
of the hexapole mode and the monopole mode. These 
symmetry-breaking comes from the rectangular grid, 
used in the FDTD computation, which does not fit in 
the hexagonal symmetry. To check this point, we re- 
calculated the far-field pattern of the hexapole mode with 
the grid of the finer grid resolution. In this computation, 
Az was fixed and only the horizontal grids (Ax, Ay) were 




Ax=Ay={1/20)(/ {1/30) (f 



(1/40) a 



FIG. 6: (a) x-polarized intensity {\Ex\'^) pattern and the 
j/-polarized intensity (\Ey\'^) pattern of the dipole mode. In- 
tensities of \Ex\'^ and \Ey\'^ are normalized by the maximum 
value of both data, (b) Calculated far-field patterns of the 
hexapole mode for the finer grid resolution {Ax, Ay). The 
vertical grid resolution Az was fixed to be o/20 throughout 
the calculation. 



changed. Fig.|Hfb) shows that, by increasing the horizon- 
tal grid resolution, the far-field pattern restores the in- 
herent six- fold symmetry. This comparison confirms that 
such unbalanced lobes are attributable to the rectangular 
FDTD grid. Note that even in case of Ax = Ay = a/20, 
all important features of the radiation can still be ob- 
tained except the slight unevenness. Thus, all the calcu- 
lations in Sec. lIVI were performed with the grid resolution 
of a/20. 



IV. ENGINEERING RADIATION PROPERTIES 
A. Defect engineering 

As discussed, the in-plane field distribution determines 
the far-field radiation pattern. In this subsection, we will 
investigate how the defect engineering method can be ap- 
plied to manipulate the radiation profile. The hexapole 
mode depicted in Sec. ITll shows null vertical emission due 
to the symmetric field distribution which has a zero dc 
component. However, as discussed in the last part of 
Sec. in it seems that the field balance can be easily bro- 
ken by introducing small structural perturbation. As 
shown in Fig. CJa), the structural tuning that we have 
applied is to enlarge two horizontal air-holes facing with 
each other. 33] From the hexapole mode field distribution 
shown in Fig. 3, it can be expected that such perturba- 
tion will break the delicate balance of the E^ distribu- 
tion. Note that noticeable amount of the dc component 
appears in the i^^-ISJI 

In Fig. Ub), the far- field patterns are plotted as we 
increase the air-hole radius. Here the same structural 
parameters used in Fig. ^ are assumed. One can clearly 
observe how the hexagonally-distributed pattern is trans- 
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FIG. 7: Radiation characteristics of the deformed hexapole 
mode, (a) Schematic of the deformed structure, where two 
horizontal air-holes facing with each other are enlarged by 
p{a). Except for this enlargement, all the structural param- 
eters are kept the same as the previous structure shown in 
Fig. 0^a). (b) Calculated far-field patterns for various values 
of the air-hole size increment p. When p = 0.04a, the far- 
field pattern shows good directionality, (c) Quality factor of 
the deformed hexapole mode as a function of the deformation 
parameter p. 



formed into the directional beaming. The radiation pat- 
tern becomes more directional with the air-hole radius. 
When the air-hole radius increment is 0.04a, the vertical 
emission looks most favorable. In this case, ~ 65% of 
the total radiated energy falls within the collection angle 
of ±30°. In other words, ^ 32.5% of the total emitted 
photons can be collected from the top. The downward 
photons could be re-directed by introducing a reflector 
underneath. In the next section, we shall discuss the ef- 
ficient unidirectional beaming strategy based on the de- 
fect tuning and the bottom reflector. The Q factors of 
the proposed hexapole mode are presented in Fig. |7Jc). 
The exponential decrease of the Q factor is observed as a 
function of the detuning parameter (p). Note that when 
the p = 0.04a, the Q factor is still not too bad, '- 20 000. 
Since the above modification of two air holes breaks 
mostly the balance of the E^ field and not that of the 
Ey field, the directional beaming obtained here will be 
almost x-polarized. Fig. |HIa) shows the polarization- 
resolved far-field patterns, where such directional beam- 
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FIG. 8; (a) a;-polarized intensity (l-B^P) pattern and the y- 



pattern of the deformed hexapole 
are normalized by the 



polarized intensity {\Ey\ ) 
mode. Intensities of \Ex\''' and \Ey 
maximum value of both data, (b) Ex, Ey, and Hz field com- 
ponents of the deformed hexapole mode whenp = 0.04a. Note 
that, although the resulting near-field distribution is almost 
identical to that of the ideal hexapole mode, the far-field pat- 
tern changes drastically. 



ing is clearly a;-polarized. With the same numerical aper- 
ture of the lens (0.5) used previously, the polarization 
extinction ratio is estimated to be ~ 70 : 1. There- 
fore, we claim that both the directionality and the linear 
polarization can be obtained from the hexapole mode 
by introducing small structural perturbation. However, 
by the same token, this result implies that the hexapole 
mode in the modified single defect cavity is not robust 
against the structural perturbation, because the electric 
fields are concentrated in the proximity of the nearest 
air-holes. This fact also explains one of the reasons why 
the observed Q factors of the hexapole mode are much 
smaller than the values obtainable in the FDTD compu- 
tation . ^ ^J Recently the other class of hexapole modes 
that are robust against the perturbation was also found, 
where the electric field energy is concentrated in the di- 
electric region between the nearest air-holes. (34 ii3(a] 

Finally, it should be noted that the near-field distri- 
bution of the deformed hexapole mode is almost identi- 
cal to that of the ideal hexapole mode and retains the 
whispering-gallery- like nature. Remember that the E^, 
Ey, and H^ field components shown in Fig. IS^b) is the 
results obtained after only very slight modification of the 
holes {p = 0.04a). It is interesting that the resulting far- 
field pattern can be changed drastically by such minute 
perturbation. Thus, as indicated by Shin et al. a thor- 
ough understanding of wavelength-size small cavities re- 
quires investigations in the both regimes; the near-field 
and the far-field. 13(1 



B. Effects of the bottom reflector 

Purcell claimed that the spontaneous emission lifetime 
is a constant that can be altered depending on the op- 
tical mode density and the electric-field intensity near 
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FIG. 9: Schematic of a novel photonic crystal unidirectional 
emitter which consists of a photonic crystal slab and a bottom 
reflector. The resonant modes are confined around the defect 
region of the slab. 



the emitter .[2i| Borrowing this concept, it can be said 
that inherent radiation properties of the PhC slab res- 
onant mode should be modified by placing a reflector 
near the cavity. ^3S| In the early days, researchers used 
PhC resonant cavities as a vehicle to study the alter- 
ation of spontaneous emission . [33, Hfl, Ell For example 
inside the periodic structure, the spontaneous emission 
can be suppressed for a certain frequency range (PBG) 
or enhanced for a particular frequency. As will be shown 
below, one can further modify the resulting radiation 
properties from the PhC cavity by using a flat-surface 
mirror. {4^ 

The PhC slab structure to be investigated here is 
shown in Fig. |51 The proposed structure is very sim- 
ple; it consists of a PhC slab and a bottom reflector. 
As a reflector, a simple dielectric substrate or a Bragg 
mirror is to be used. Such a structure is realistic since 
the thickness of each layer can be controlled during the 
wafer growth. We note that similar planar systems were 
studied by Schubert[43| and Benisty 44] in an attempt 
to alleviate the poor light extraction efficiency of micro- 
cavity LEDs. To investigate effects of the presence of 
the bottom reflector, we have calculated the change of 
Q factors and far-field radiation patterns as we vary the 
air-gap distance between the PhC slab and the bottom 
reflector (see Fig. llO|) . Here, the Q factor represents the 
lifetime of the resonant mode (r = Qjoj). As a PhC 
resonant mode, the quadrupole mode shown in Fig. |Slis 
used. As a bottom reflector, we assumed a simple dielec- 
tric substrate whose refractive index is 3.4. 

Compared with the original far-fleld pattern shown 
previously in Fig. [S] the radiation proflle shown in 
Fig. EH changes noticeably with the air-gap size. The 
far-field patterns show the signature of concentric ring- 
shape modulation laid upon the original far-fleld pattern, 
which is reminiscent of the Fabry-Perot fringes 45] in 
the Fabry-Perot etalon. Although the reflectivity of the 
bottom GaAs (or InP) substrate is only '^ 30%, the in- 



terference between the original wave and the reflected 
wave is easily noticeable. Upon inspection of the far- 
field patterns in Fig. llOf a). one can find that the sim- 
ilar far-fleld patterns show up again after the air-gap 
increment of ^ 1.75a (for example, compare 1.75a and 
3.50a). This increment corresponds to around a half- 
wavelength (~ 1.75a X 0.3179A/a ~ 0.55A). It should be 
noted that the directional beam cannot be achieved from 
the quadrupole mode. To achieve unidirectional emis- 
sion, we need to investigate the inherently directional 
resonant mode such as the deformed hexapole mode. 

Similar systematic dependence on the air-gap size is 
also observed in the Q factor of the resonator, at a period 
of approximately A/ 2. Note that the Q factor decreases 
drastically when the air-gap size becomes less than 0.25A. 
There are two distinct reasons for this phenomenon; one 
is the TE-TM couphng loss through the PhC slab H^ 
and the other is the diffraction loss in the downward di- 
rection. The TE-TM coupling loss always occurs when 
the vertical symmetry of the PhC slab structure is bro- 
ken. The reason for the increase of the diffraction loss 
to the downward direction is that the light-cone of the 
bottom substrate is effectively enlarged by a factor of n 
(the refractive index of the substrate). When the bot- 
tom substrate is too close from the slab, the tail of the 
evanescent fleld begins to feel the substrate. Such an 
expanded light-cone will allow the more plane-waves to 
be coupled to substrate-propagating modes. Therefore, 
in order to prevent this loss, the air-gap size larger than 
A/2 is generally recommended. It is interesting to find 
that the Q factor can be enhanced even in the presence 
of the bottom substrate. In fact, when the air-gap size 
is larger than A/2, the Q factor oscillates around the 
original value. Note that the interference phenomena ob- 
served here are similar to the cases of the dipole antenna 
placed near an ideal plane mirror. 42] 

The case of the deformed hexapole mode is worth spe- 
cial attention. Fig. Illf a) shows the far-fleld radiation 
patterns with a perturbation parameter p = 0.04a (see 
Fig. [7||. The same bottom substrate of index 3.4 is used. 
Here again, one can see that the original far-field pat- 
terns are modified by concentric ring-shape modulations. 
Comparing the patterns of 1.25a and 3.00a, one can wit- 
ness that the similar vertical enhancement is restored 
after the air-gap is increased by additional thickness of 
~ 1.75a. This air-gap thickness corresponds to, approxi- 
mately, a half-wavelength (~ 1.75a x0.3192A/a ^ 0.56A). 
At this point, the vertical enhancement condition seems 
to be summarized as 2d ~ mX , where d is the air-gap 
size and m integral multiple. Moreover, the air-gap size 
of 3.00a where the good directionality is obtained is ap- 
proximately equal to one wavelength {^ 0.96A). If one 
wants to find the rigorous conditions of optimal verti- 
cal enhancement, the analyses need to include the effects 
of the optical thickness of the PhC slab and the phase 
change upon reflection. In the next section, we shall de- 
velop a simple, plane-wave based interference model to 
explain the vertical enhancement conditions. From the 
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FIG. 10; Radiation characteristics of the quadrupole mode 
in the presence of a simple dielectric substrate with a refrac- 
tive index of 3.4. (a) Calculated far-field patterns for various 
values of the air-gap thickness. Here, each intensity value has 
been normalized to the maximum value of its own data, (b) 
Calculated quality factors as a function of the air-gap size. A 
horizontal dotted line represents the original quality factor in 
the absence of the reflector. 
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FIG. 11: Radiation characteristics of the deformed hexapole 
mode in the presence of a simple dielectric substrate with 
a refractive index of 3.4. (a) Calculated far-field patterns 
for various values of the air-gap thickness. Here, each inten- 
sity value has been normalized to the maximum value of its 
own data, (b) Calculated quality factors as a function of the 
air-gap size. A horizontal dotted line represents the original 
quality factor in the absence of the reflector. 



Q factors, one can clearly see, again, the modification 
of the intrinsic radiation lifetime, where the oscillation 
period around the dotted line is ^ A/2 [See Fig, illf b^]. 
Note that there seems to be no definite relation between 
the Q factor optimization and the vertical enhancement 
condition; The optical loss of the cavity (~ l/Q) is a re- 
sult of the sum of the emitted power {J{dP/dfi)di^). A 
simple guess of the Q by the emission pattern itself can 
be misleading. 



In principle, one can completely block the downward 
propagating energy by using a highly-reflective Bragg 
mirror. For example, a stack of 25 AlAs/GaAs pairs 
give a reflectivity in excess of 0.998. [23 In the subsequent 
analyses, we will use a Bragg mirror made of a stack of 20 
AlAs/GaAs pairs whose reflectivity is ^ 99%. Remem- 
ber that the FDTD grid resolution for the z direction 
[(l/20)a] is, in fact, not sufficiently fine to fully describe 
the each Bragg layer (around 10 grids are used to describe 
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FIG. 12: Radiation characteristics of the deformed hexapole 
mode in the presence of a Bragg mirror, (a) Calculated far- 
field patterns for various values of the air-gap thickness. Here, 
each intensity value has been normalized to the maximum 
value of its own data, (b) Calculated quality factors as a 
function of the air-gap size. A horizontal dotted line repre- 
sents the original quality factor in the absence of the reflector. 



time by 75% in comparison with that without a reflector. 
Note that, except for such increased interference effects, 
all other characteristic features in the far-field patterns 
still remain almost unchanged. In short, we confirmed 
that both the radiation lifetime and the radiation pat- 
tern can be controlled by simply adjusting the distance 
between the PhC slab and the reflector. As shown in 
the far-field pattern of 3.00a, one can achieve a sharp 
directional emission in which most of the emitted power 
(^ 80%) is concentrated within ±30°. 



C. Design rule for directional emission 

In the preceding subsection, we have shown that, by 
using the 3-D FDTD method and the near- to far-field 
transformation algorithm, the radiation patterns of the 
PhC defect modes can be altered by the bottom reflector. 
Especially the choice of ^^ lA air-gap seems to give the 
best directional emission. In this subsection, we will try 
to find conditions for the directionality, by developing a 
simple analytic model based on plane- wave interferences. 

We assume that the PhC slab can be approximated 
as a uniform dielectric slab having an effective refractive 
index[47j and the PhC resonant mode can be treated as 
a point dipole source embedded in the middle of the slab 
and the radiations are described by scalar plane- waves. 
We also assume that the bottom reflector is positioned 
sufficiently far (> A/2) from the PhC slab. The reflector 
does not perturb too much the original characteristics of 
the resonant mode. Thus, the initially upward and down- 
ward propagating waves (See Fig. lCT a)) can be assumed 
to be symmetric, 1 and 1. The downward propagating 
component will eventually be redirected upward to make 
an interference with the upward propagating component. 
Although these assumptions seem to be very simple, sim- 
ilar approaches based on the plane-wave interferences 
have already been used for the calculation of radiation 
patterns of atomic dipoles in a planar microcavitv. ^44. 48l | 
As indicated by Dowling4Sj, the mode structure of the 
electromagnetic fleld is a classical phenomenon. 

In the simplest case when multiple reflections between 
the PhC slab and the bottom reflector can be neglected 
|49| . a constructive inteference condition at the vertical 
becomes 



one pair of AlAs/GaAs).i2ll| As a result, the reflectivity 
would be slightly smaller than the actual value. 

Far-field patterns in presence of a Bragg mirror are 
collectively shown in Fig. I12f a'l. In comparison with 
Fig. IllT a). the concentric ring-shape modulation is more 
pronounced due to the high reflectivity of the bottom 
structure. In the Q factor plot shown in Fig. I12f b). the 
modification of the radiation lifetime is also prominent. 
The maximum Q factor is increased from ~25 000 to ^35 
000, which represents an increase of the radiation life- 



^2z<pgzeg»0 ^ ^ 



(18) 



where 2(/p is a round-trip phase in the air-gap. e and 
are the phase changes at the bottom reflector |5(]l | 
and through the PhC slab, respectively (See Fig.^fa)). 
When the radiation into an angle 9 is considered. 



Lp — [2ix IX)dcos6 



^ = (27rne///A)L^l-(l/n2^^)si 



sin^6' 



(19) 
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FIG. 13: Vertical emission in the presence of the perfectly conducting mirror (PCM), (a) Simple vertical interference model. 
The emission from the cavity is divided by the upward and the downward propagating waves, (b) Calculated vertical emission 
intensity \Sf {6 = 0) from the initially downward propagating wave. Here, € — tt and tp = 2tv are assumed, (c) Calculated 
vertical emission enhancement factor W = |1 + S\'^ {6 = 0) as a function of the slab thickness for various values of the air-gap 
thickness ranging from l.OOA to 0.85A. 



where d and L are the air-gap thickness and the slab 
thickness, respectively. Considering the typical thickness 
of the PhC slab of « A/(2n), one may assume (/) « tt.^I 
Then, Eq. l(TF|) becomes 



exp(— i47r(i/A) « 1 



(21) 



This condition is consistent with the previous obser- 
vation {2d ~ mX) which has been deduced from the 
far-field patterns. However, in general, the PhC slab 
thickness does not satisfy this 'slab resonance condition 
(L = A/(2n))' and the residual reflection can not be com- 
pletely cancelled. Then, what should be changed in the 
vertical enhancement conditions [Eq. (|18|l ] when multiple 
reflections become significant? 

Fig- EI a-) describes the situation where the emission 
from the cavity is divided by the upward and the down- 
ward propagating waves. Here, we focus only on the ver- 
tical interference conditions {9 = 0). ro and to are coeffi- 
cients of amplitude reflection and transmission for a sin- 
gle dielectric interface, respectively. Let S be the sum of 



all the waves detected upward which initially propagate 
downward. Then, an interference between the originally 
upward propagating wave and the S will be described by 
1 + S", and the resultant intensity will be proportional to 
|1 -|- S'p. To sum all the waves which undergo infinitely 
many multiple reflections in the slab and in the air-gap, 
we introduce a simple method which take advantage of 
the similarity of such multiple reflections and transmis- 
sion. 

We trace several reflections and transmissions of the 
wave which start out at point A in the initial stage. At 
point B, the wave is divided into two. Here, the down- 
ward propagating component will undergo the same mul- 
tiple reflections as the wave at A. Let Q he the sum of 
all the waves detected upward which propagate upward 
at B, then we get the following equation for the S. 



S = e^'^^e' {toQ + {-ro)S) 



(22) 



Now, let us proceed the upward propagating wave at B 
more. The wave once reflected at the top of the slab will 
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condition for a constructive interference becomes 



FIG. 14: Calculated far-field radiation of the deformed 
hexapole mode in the 'A' air-gap condition, (a) Far-field pat- 
terns for various values of the slab thickness ranging from 
0.6a to 0.9a, which approximately correspond to from ~ vr to 
~ 1.27r. (b) Polar plot of the far-field pattern. Note that the 
far-field pattern becomes more directional as the slab thick- 
ness increases. From the far-field pattern of L = 0.9a, it is 
found that more than 80% of total emitted power is directed 
within a small divergence angle of ±30°. 



reach at point C, there, two divided waves have resem- 
blance to the waves which contribute to the Q and the 
S, respectively. Describing these processes, 



Q = e'"^ (to + roe'^roQ + to^)) (23) 

We get the S by equating Eq. ^ and Eq. f^ . 



j-2 # 



(1 - r^e2''^)(ro -l- g-^'Ve-^e) - rotle'^'-'l' 



(24) 



Now, let us examine two representative cases of the 
PhC slab; 'slab resonance' and 'slab antiresonance'. 



TT, I.e., e"^ — —1. From the well 



(1) Slab Resonance 

This is a case when (j) = tt, i.e. , e*"^ 
known result for the transmittance in a uniform dielectric 
slab, 1231 we know that the slab becomes transparent and 
multiple reflections in the air-gap do not occur. Thus, we 
can easily guess that this case restores the simplest two- 
beam interference result dealt in Eq. H18|l . By inserting 



oi0 



-1 into the Eq. if^ . we get 



S = 



ig(-l) 



-1 



g-2^■^g-^£ 



g — 2z</Pg — ie 

Resultant radiation intensity becomes 



W= ll + S'P 



1 



g-2»Vg-*£|2 



(25) 



(26) 



This is the famous two-beam interference formula, if 
the phase change at the bottom mirror is e = tt, then the 



2(p = 2rmr 



The maximum radiation intensity is W, 
the above condition. 



(27) 

max = 4 with 



(2) Slab Antiresonance 
This is a case when 



1.57r, i.e., e 



i4> 



The 



reflectance from the dielectric slab is maximized with this 



S ■ 



Inserting e"^ 


-itl 


into Eq. ^ 


(l + r2)(ro + 


g-2z^ 


e-'') + rotl 



(28) 



Here, we let e again be tt, then one can easily prove that 
the condition of the air-gap to maximize \S\ is e"^"'' = 1 
(air-gap resonance). 



Dr, 



-iti 



(1 



rDira 



1) 



(1 



rotl 



.(l + ^o) 



(l-ro)2 (1-ro) 



(29) 



Assuming the effective refractive index of the slab 
(rieff) by 2.6, then tq = {ueff - l)/{neff + 1) « 0.44. 
The vertical radiation intensity becomes 



W=\l + S\^ = v 



l + rp 
l-ro 



7.9 



(30) 



which can be larger than that of the slab resonance. 
Note that the condition for maximizing W is slightly 
different from the slab antiresonance, because of a phase 
difference between two complex numbers {Smax is a pure 
imaginary while '1' is a real). 



In Fig. liaf c'). we plot W a,s a. function of the slab thick- 
ness {(p) for various air-gap sizes ((/?). A blue line repre- 
sents the result when the air-gap size exactly equals the 
emission wavelength {d — l.OOA). Other lines represent 
the air-gaps as indicated. Here, it is found that 1) the 
overall emission intensity decreases as the air-gap is de- 
tuned from the 'A' condition and 2) when the air-gap size 
is 'A'(on-resonance), there exists a relatively wide con- 
cave plateau in the neighborhood of the 'slab resonance' 
(0 = tt). In this case, the vertical emission is enhanced 
by a factor larger than 4, over a wide range of slab thick- 
nesses. Although the better vertical enhancement can be 
obtained near the 'slab antiresonance' {(j) = I.Stt), the 
thickness of the PhC slab is better to be chosen near the 
'resonance' condition, because the thicker slab tends to 
support unwanted higher-order guided modes., 2t|| Thus, 
thanks to the multiple reflection effect, the vertical en- 
hancement condition now becomes rather simple. Above 
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FIG. 15: Electric-field intensity distributions of the deformed 
hexapole mode detected in the x-z plane (a) in the absence 
of the reflector (b) in the presence of the Bragg mirror. 
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simply adjusting the air-gap size to be A, one can obtain 
fairly good directional patterns. The directionality 
improves continuously until the slab thickness becomes 
near the slab antiresonance. In a practical point of 
view, the air-gap size is a fixed quantity which cannot 
be changed after the wafer growth [53|. Thus, the most 
important design rule for the directional emission is to 
choose the air-gap size (thickness of a sacrificial layer) to 
be a target wavelength 'A'. Then, the remaining critical 
issue will be that of the lithographic tuning that matches 
the resonant wavelength of the mode of interest with the 
emission wavelength of quantum well or quantum dots. 



V. SUMMARY 



analyses suggest that the focus should be placed more 
on the air-gap size that has be exactly matched with the 
emission wavelength A than on the PhC slab thickness 
that is already designed near the antiresonance condi- 
tion. 

To confirm the above argument, we have calculated 
far-field patterns of the deformed hexapole mode using 
the 'A' air-gap. We have increased the slab thickness 
from 0.6a to 0.9a, which correspond to approximately 
the phase thickness from ~ tt to ~ 1.27r [SJ. Since reso- 
nant mode frequencies are changed by the slab thickness, 
the air-gap sizes have been adjusted to fit the 'A' con- 
dition. As shown in Fig. ^J good directional patterns 
are obtained. It is interesting to note that the far-field 
pattern becomes more directional as the slab thickness 
increases. Therefore, as expected in the previous argu- 
ment, the slab slightly thicker than that satisfying the 
resonance condition is preferable for the enhanced ver- 
tical beaming. From the far-field pattern of L = 0.9a, 
it is found that more than 80% of the emitted power is 
contained within an angle of ±30°. 

Figure ^| shows the electric-field intensity distribu- 
tions detected in the x-z plane. For comparison, the 
case of a symmetric PhC slab cavity without a bottom 
reflector is also shown. In Fig. llSf b). one can clearly 
see that the downward propagating waves are effectively 
redirected by the Bragg mirror and the unidirectional 
emission is achieved. In other words, we have shown that 
inherent radiation patterns of the PhC slab resonant 
mode can be controlled either by modifying the defect 
structures, or by placing a reflector near the cavity. By 



We have shown that the radiation pattern and the ra- 
diation lifetime of the 2-D PhC slab resonant mode can 
be controlled by placing a mirror near the cavity. To cal- 
culate the far-field pattern, we have developed a simple 
and efficient simulation tool based on the FDTD method 
and the near- to far-field transformation. We have dis- 
cussed two distinct physical mechanisms responsible for 
the modification of the far-field radiation pattern. Ei- 
ther by tuning the defect structure or by placing a re- 
fiector near the cavity, we have shown that the unidi- 
rectional emitter whose energy is concentrated within a 
small divergence angle (±30°) can be achieved. Simple 
plane-wave based interference model has revealed that 
the vertical enhancement condition is satisfied when the 
air-gap size is equal to the emission wavelength. To fur- 
ther enhance the directionality, the slightly thicker PhC 
slab design than the 'slab resonance' condition would be 
recommended. The proposed structure is practical since 
the critical design parameter can be precisely controlled 
during the epitaxial growth. 



Acknowledgments 

One of the authors, S. H. Kim, would like to thank S. 
H. Kwon, M. K. Seo, J. K. Yang for their valuable assis- 
tance and discussions. S. H. Kim would also like to thank 
J. M. Gerard in CEA in France for his valuable discus- 
sions. This research was supported by a Grant from the 
Ministry of Science and Technology (MOST) of Korea. 



[1] J. D. Joannopoulos, P. R. Villeneuve, and S. Fan, Nature 
(London) 386, 143 (1997). 

[2] S. G. Johnson, S. Fan, P. R. Villeneuve, J. D. Joannopou- 
los, and L. A. Kolodziejski, Phys. Rev. B 60, 5751 (1999). 

[3] S. H. Kim, G. K. Kim, S. K. Kim, H. G. Park, Y. H. Lee, 
and S. B. Kim, J. Appl. Phys. 95, 411 (2004). 

[4] O. Painter, R. K. Lee, A. Yariv, A. Scherer, J. D. 



O'Brien, P. D. Dapkus, and I. Kim, Science 284, 1819 

(1999). 
[5] H. G. Park, J. K. Hwang, J. Huh, H. Y. Ryu, Y. H. Lee, 

and J. S. Kim, Appl. Phys. Lett. 79, 3032 (2001). 
[6] H. Y. Ryu, S. H. Kwon, Y. J. Lee, Y. H. Lee, and J. S. 

Kim, Appl. Phys. Lett. 80, 3476 (2002). 
[7] Y. Akahane, T. Asano, B. S. Song, and S. Noda, Nature 



14 



(London) 425, 944 (2003). 
[8] H. Y. Ryu, M. Notomi, and Y. H. Lee, AppL Phys. Lett. 

83, 4294 (2003). [31] 

[9] T. Yoshie, A. Scherer, J. Hendrickson, G. Khitrova, H. M. 

Gibbs, G. Rupper, C. Ell, O. B. Shchekin, and D. G. [32] 

Deppe, Nature (London) 432, 200 (2004). [33] 

[10] J. M. Gerard and B. Gayral, Proc. SPIE 5361, 88 (2004). [34] 

[11] J. Vukovic and Y. Yamamoto, Appl. Phys. Lett. 82, 2374 [35] 

(2003). [36] 

[12] K. Hennessy, A. Badolato, P. M. Petroff, and E. Hu, 

Photonics and Nanostructures (Elsevier) 2, 65 (2004). [37] 

[13] H. G. Park, S. H. Kim, S. H. Kwon, Y. G. Ju, J. K. Yang, [38] 

J. H. Baek, S. B. Kim, and Y. H. Lee, Science 305, 1444 

(2004). [39] 

[14] S. Noda, M. Imada, M. Okano, S. Ogawa, M. Mochizuki, 

and A. Chutinan, IEEE J. Quantum Electron. 38, 726 [40] 

(2002). 
[15] T. Asano, M. Mochizuki, S. Noda, M. Okano, and [41] 

M. Imada, IEEE J. Lightwave Technol. 21, 1370 (2003). [42] 

[16] S. Noda, A. Chutinan, and M. Imada, Nature (London) 

407, 608 (2000). 
[17] H. Takano, Y. Akahane, T. Asano, and S. Noda, Appl. [43] 

Phys. Lett. 84, 2226 (2004). 
[18] H. G. Park, J. K. Hwang, J. Huh, H. Y. Ryu, S. H. Kim, [44] 

J. S. Kim, and Y. H. Lee, IEEE J. Quantum Electron. 

38, 1353 (2002). [45] 
[19] H. Y. Ryu, H. G. Park, and Y. H. Lee, J. Sel. Top. Quan- 
tum Electron. 8, 891 (2002). [46] 

[20] L. A. Coldren and S. W. Corzine, Diode Lasers and Pho- 
tonic Integrated Circuits (Wiley, New York, 1995). [47] 

[21] A. Taflove and S. C. Hagness, Computational Elec- 
trodynamics: the finite- difference time-domain method [48] 
(Artech House, Norwood, MA, 2000), 2nd ed. 

[22] J. Vukovic, M. Loncar, H. Mabuchi, and A. Scherer, [49] 
IEEE J. Quantum Electron. 38, 850 (2002). 

[23] H. Y. Ryu, J. K. Hwang, and Y. H. Lee, J. Appl. Phys. 

88, 4941 (2000). [50] 

[24] O. J. Painter, A. Husain, A. Scherer, J. D. O'Brien, 
I. Kim, and P. D. Dapkus, IEEE J. Lightwave Technol. 
17, 2082 (1999). 

[25] K. Srinivasan and O. Painter, Opt. Express 10, 670 [51] 
(2002). 

[26] S. H. Kim and Y. H. Lee, IEEE J. Quantum Electron. 

39, 1081 (2003). 

[27] K. Sakoda, Optical Properties of Photonic Crystals [52] 

(Springer- Verlag, Berlin Heidelberg, 2001). 
[28] J. D. Jackson, Classical Electrodynamics (Wiley, New 

York, 1974). 
[29] W. H. Press, B. P. Flannery, S. A. Teukolsky, and W. T. 

Vetterling, Numerical Recipes in C: The Art of Scientific 

Computing (Cambridge University Press, 1992), 2nd ed. 
[30] D. J. Shin, S. H. Kim, J. K. Hwang, H. Y. Ryu, H. G. 



Park, D. S. Song, and Y. H. Lee, IEEE J. Quantum Elec- 
tron. 38, 857 (2002). 

M. L. Povinelli, S. G. Johnson, J. D. Joannopoulos, and 
J. B. Pendry, Appl. Phys. Lett. 82, 1069 (2003). 
O. Painter and K. Srinivasan, Opt. Lett. 27, 339 (2002). 
S. K. Kim, (unpublished). 
S. H. Kwon, (unpublished). 
G. H. Kim, (private communication). 
S. K. Kim, J. H. Lee, S. H. Kim, I. K. Hwang, Y. H. Lee, 
and S. B. Kim, Appl. Phys. Lett. 86, 031101 (2005). 
E. M. Purcell, Phys. Rev. 69, 681 (1946). 
G. Bjork, S. Machida, Y. Yamamoto, and K. Igeta, Phys. 
Rev. A 44, 669 (1991). 

J. K. Hwang, H. Y. Ryu, and Y. H. Lee, Phys. Rev. B 
60, 4688 (1999). 

R. K. Lee, Y. Xu, and A. Yariv, J. Opt. Soc. Am. B 17, 
1438 (2000). 

H. Y. Ryu and M. Notomi, Opt. Lett. 28, 2390 (2000). 
E. A. Hinds, tn Cavity Quantum Electrodynamics (Aca- 
demic Press, Inc, Orlando, 1994), edited by P. R. 
Herman. 

E. F. Schubert, Y. H. Wang, A. Y. Cho, L. W. Tu, and 
G. J. Zydzik, Appl. Phys. Lett. 60, 921 (1992). 
H. Benisty, H. D. Neve, and C. Weisbuch, IEEE J. Quan- 
tum Electron. 34, 1612 (19982). 

M. Born and E. Wolf, Principles of Optics (Cambridge 
University Press, 1980), pp. 329,330. 
Y. Tanaka, T. Asano, Y. Akahane, B. S. Song, and 
S. Noda, Appl. Phys. Lett. 82, 1661 (2003). 
S. Fan and J. D. Joannopoulos, Phys. Rev. B 65, 235112 
(2002). 

J. P. Bowling, M. O. Scully, and F. DeMartini, Opt. 
Commun. 82, 415 (1991). 

When an optical thickness of a PhC slab is 'slab res- 
onance condition', the PhC slab becomes transparent 
(zero reflectance). 

When we use a Bragg mirror as a bottom reflector, a 
reflection phase is n at the Bragg frequency. This phase 
shift is the same with the case of a perfect conducting 
mirror. 

Optical thicknesses of the PhC slab are estimated from 
the transmittance spectrum which has been obtained by 
using the 3-D FDTD method with periodic boundary 
conditions. For references, see Fan et al. |4^ 
Recently, LEOM group in France has demonstrated PC- 
MOEMS, where one can control the air-gap size by using 
electro-static force. See, for instance, J. L. Leclercq, B. 
BenBakir, H. Hattori, X. Letartre, P. Regreny, P. Rojo- 
Romeo, C. Seassal, P. Viktorovitch, Proc. SPIE, 5450, 
300 (2004). 



